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Abstract. We present a model of an ATP-fueled molecular machine which push 
a polymer through a pore channel. The machine acts between two levels (working- 
. waiting), and the working one remains active for a hxed time giving a constant 

£Zj ' force. The machine activation rate can be put in relationship with the available 

ATP concentration in the solution, which gives the necessary energy supply. The 
translocation time shows a monotonic behaviour as a function of the activation 
frequency and the velocity follows a Michaelis-Menten law that arises naturally in 
this description. The estimation of the stall force of the motor follows a corrected 
| Michaelis-Menten law which still is to be checked in experimental investigation. The 

results presented agree with recent biological experimental findings. 
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1. Introduction 



In the last years there has been an important progress in the experimental and 
theoretical study of transport mechanisms of molecules inside cells and/or through 
cell membranes [1]. This effort is attracting nowadays more and more attention from 
researchers in different scientific disciplines. The increasing interest in this subject 
is related to the intrinsic importance of understanding the basic mechanisms of the 
living systems, but also to the enormous improvement of technological capabilities 
at nanometer length scale. These, on the one hand, allow to detect and to measure 
mechanisms at the nanoscale, and, on the other hand, open the possibility of 
constructing from scratch structures (with both natural and synthetic materials) able 
to imitate the biological functioning [2j. In this context, the passage of biomolecules 
through nanopores is ubiquitous, both in biological and nanotechnological processes. 
Examples of these two types are the passage of mRNA through nuclear pores j3] or the 
translocation of DNA in graphene pores [I] . 

In most cases, translocation is driven by constant fields in the pore or by the 
difference in chemical potential between the two sides of the membrane. However, in 
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some cases the translocation is assisted by an ATP-based molecular motor [5]. This is 
the case of DNA bacteriophages in which the incoming DNA has to overcome a huge 
pressure inside the virus capsid. This makes this kind of motors, possibly, the most 
powerful of those known. In this type of motors ATP hydrolysis is the fuel of the 
process. 

In this paper we model the translocation process of a Id chain pushed by a molecular 
motor activated by ATP absorption. The motor is able to drive with a constant force a 
polymer chain in one direction, while in its activated state. The polymer diffuses freely 
otherwise [61 [7] . The work reveals the Michaelis-Menten (MM) behaviour of the polymer 
velocity and relates it with the MM enzymatic reaction, according to a microscopic re- 
interpretation of the MM kinetics which is a very actual topic of investigation [HI [91 [10] . 

We will also study the behaviour of the motor against a pulling force, the motor 
stall force, and its ability to package the polymer in a finite region, as a first approach 
to a capsid effect. 

The aim of this work is to present a simple model which captures the main 
physical ingredients of the process. Remarkably, the model is able to well describe 
some experimentally observed results and make new predictions. The model could be 
also applied to different kinds of systems. In this spirit, we do not pretend here to give 
a detailed description of a particular molecular motor. 



2. The model 



The polymer is modeled as a 1-dimensional chain of N dimensionless monomers 
connected by harmonic springs jllj. The total potential energy is Vhar = § Sili^i+i ~~ 
%i — do) 2 , where k is the elastic constant, Xi the position of the i-th particle, and do the 
equilibrium distance between adjacent monomers. Along the work do = 1 and k = 1. 
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Figure 1. Scheme of the dichotomous pushing force acting on a polymer chain 
compound by N monomers. Tm is the working time of the motor, supposed fixed. 
T is the mean waiting time in the inactive state. The motor acts on a region of length 
Lm, where the motor action can be expressed by the linear potential —i](t)FMX (with 
i] = 1 for the active state and r\ = for the inactive one), indicated in the middle 
figure by two dotted lines. 
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The translocation is helped by the presence of an ATP activated molecular motor. 
It has spatial width L M and we set at x = its right edge (see figured]). The motor is 
characterized by a fixed working time Tm which follows the ATP absorption, which, in its 
turn, occurs after a mean waiting time T, depending on the very ATP concentration. In 
this sense the ATP molecules act on the motor as a shot noise contribution able to switch 
on its activity. The motor exerts a dichotomous force F^^(t) on the particles inside 
the motor (x G [— Lm,0]), where rj(t) is 1 during the working time and otherwise. As 
shown in [12], the ATP absorption follows an exponential distribution of waiting times: 
the probability for an ATP adsorption in a time between t' and t'+dt' after the last 
activity is proportional to e~ t '^ T . Thus, the activation probability of the motor is given 
by P t i — 1 — e~ l> l T with t! the motor residence time in its inactive state. 

The dynamics of the i th monomer of the chain is described by the following 
overdamped Langevin equation: 

Xi = -V^ r + F MVi (t)+^(t), (1) 

where the damping (considered the same for all the monomers, which feel the viscosity 
independently on each other) is included in the normalized time. In the above 
dimensionless equation d Q is taken as the unit of length, the force is expressed in 
units of kdo, and the energy in units of kd 2 Q . are the Gaussian uncorrected 

thermal fluctuations which follow the usual statistical properties (&(£)) = and 
(Ut)Ut + t)) = 2D5 tJ 5(r) with = 1...N). 



3. Energy supply 

The energy used by the motor to provide the driving is given by the ATP hydrolysis. In 
our model, each ATP hydrolysis just activates the motor for a fixed working time Tm- 
This description avoids more complex details which are beyond the spirit of our work. 

As a consequence of the ATP absorption, the motor changes its conformational state 
and exerts a mean force Fm during a fixed time Tm = l/^o- After that time the motor 
returns to its inactive state. A new force is applied when the next ATP suitable quantity 
is absorbed. It happens after a mean time T = 1/u, which follows an exponential 
distribution of waiting times. That way in our description the activation frequency 
is proportional to the ATP concentration [13]: v oc [ATP]. The hypothesis that the 
motor works for a fixed time and that the statistics of the arrival of the ATP molecules 
happens in an exponential distribution appears realistic in good approximation as put 
in evidence by different experimental works P [101 [El El [15] . 

With this definition the Michaelis-Menten law arises naturally from the model. The 
motor duty ratio, the fraction of time that it is active, is given by 
Tm v [ATP] 

T M + T v + is k M + [ATPy [) 
which is a Michaelis-Menten law. The only hypothesis included in the derivation of the 
last equality is that the motor activation rate v is proportional to the ATP concentration, 
being kM — [ATP]i> Q /i> the Michaelis constant. 
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The relationship between the mechanical description presented and the Michaelis- 
Menten law is deeper than only the statistical ansatz v oc [ATP]. In the Michaelis- 
Menten enzymatic reaction, 



the rate k\ represents the probability to form the compound Z per unit of time and 
per unit of S ([ATP]), and k<i gives the probability to form the product P per unit 
of time. In our mechanical and individual case (single motor and single ATP event) Z 
represents the ATP bound to the motor, which occurs with the frequency v (~ k\ [ATP]), 
while P represents the motor action, which is completed within a time Tm — V^o 
(z^o ~ ^2). These relationships are in agreement with the definition of the Michaelis 
constant ku = k 2 /ki. 

4. Results 

We have numerically solved the Langevin equation ([Q). We are mainly interested in the 
behaviour of the translocation time r and velocity v as a function of the motor activation 
frequency v. We will compare our results to some theoretical predictions, derived below. 
For every computed point we performed at least N exp = 10, 000 numerical experiments 
(for low frequency up to 50, 000) using a stochastic Runge-Kutta algorithm with a time 
step dt = 0.01. The polymer is compound by N monomers and starts with all the spring 
at the rest length and the last monomer of the chain at xn = 0, just at the exit of the 
motor. The noise intensity is held fixed at the value D = 0.01 and the intensity of the 
force is F M = 0.1. We have chosen L M /d = 5.5. 

Polymer translocation. We start studying the mean translocation time and velocity 
of the polymer driven by the motor. With respect to the polymer velocity, by summing 
up the terms of equation (CQ) and averaging in time, we obtain for the centre of mass 
of the chain the mean velocity v: 



Here n™ ot (v) is the mean number of monomers inside the motor during its activity, 
a number which is expected to depend weakly on v. Thus, the mean velocity of the 
polymer depends on the force felt by the n^ ot monomers inside the machine which 
operates for fraction of time vj (vq + u), and the polymer shows a Michaelis-Menten law 
for the velocity weakly moderated by the function n™ ot (i > ). Note that this velocity goes 
to zero as 1/N for a large chain as expected for a motor which acts on a small number 
of monomers and a polymer which moves in a dissipative media. 

We also want to mention that equation it is also valid in other interesting 
situations too. For instance, systems with random distributions of Tm values with small 
dispersion around its mean value or the case of a minimum threshold time for the motor 
activation, which effectively imposes a maximum cutoff frequency in the system (in this 
case 1/v = T thres + T). 




(3) 




(4) 




Figure 2. Log-Log plot of the mean first passage time r as a function of the mean 
frequency of the fluctuating force. Upper inset: mean velocity, following the Michaelis 
Menten law given by equation ([5]). Bottom inset: mean number of particle inside the 
motor during its active state. The thermal noise intensity is D = 0.01 and k = 1. 



Figure [2] shows the main observables of the system and their frequency dependence. 
The translocation time r is computed as the mean first passage time (MFPT), i.e. the 
average over the N exp realizations of the time spent by the centre of mass of the chain 
to reach the position x = 0. It is observed that this time decreases monotonously as v 
increases and reaches a limit value for large enough values of v (T goes to zero and the 
motor is 'on' most of the time). The relation between translocation velocity v and time 
t is not trivial at all, as can be seen in [16] and |17j . 

The physics of the problem is regulated by the mean number of monomers in the 
motor during the working time, n^ ot , a number directly related with the elasticity of the 
chain. As seen in the inset, for the used parameters this number is close to L M /d and 
does not change importantly in all the frequency range covered. Thus, as predicted by 
equation (TjJ the polymer velocity follows a Michaelis-Menten law (see inset in figure |2]): 

v * v° HR /(l + is M (5) 

with v RR the high frequency limit (1^^=0.0450 in the plotted case, which is close to 
0.0458=F M n M /iV with F M = 0.1, N = 12 and n M = L M /d = 5.5). Equation © allows 
for a direct experimental fit once the velocity at high ATP concentration (the high rate 
limit of our system) is measured. With respect to the low frequency limit, we can see 
from the equation that v goes to zero as v ~ v° hr {u/uq) ~ {F m um / N){v / u Q ) . 

Elastic Constant. In order to better understand how the elasticity of the chain 
acts on the translocation dynamics, we have computed the mean velocity and first 
passage time (v and r) for different values of the elastic constant between monomers 
k (see figure [3]). We notice that similar behaviour is observed for the different values 
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Figure 3. Translocation time and velocity (inset) for different values of the elasticity 
of the chain k. Lines in the inset show the fit to v = Vjjh/(1 + bv$/v) where b = 1 for 
k > 0.5. 



considered. As expected, a smaller k gives a lower velocity, the polymer enlarges and 
the mean number of monomers in the machine diminishes. At low k, n™ ot changes 
significantly with v and v shows deviations with respect to the predictions of equation (jSJ) 
but can be nicely fit by a slightly modified law (see caption of figure [3]). For high k 
(k > 1) we reach the rigid chain limit of the system. 

Translocation in the presence of a pull force. A set of simulations has been 
performed by applying a pull force F p on the left extremum of the chain, acting against 
the motor (see figure H]). The initial condition is set with the polymer centre of mass in 
the centre of the machine. The velocity of the centre of mass is measured waiting for 
the exit on the left or on the right of the potential region. 

In this case, the mean velocity is given by 
= n™ ot (is; F p ) _ Fp 

N 1 + uo/u N ' 1 ' 

where the last term F p /N modifies equation (j4]) by taking into account the presence of 
an external force acting against the machine. There we emphasize that the presence of 
a pulling force also modifies the value of n£f ot with respect to the unforced value. 
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Figure 4. A pull force F p is applied at the chain to measure the stall force of the 
motor. 
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Figure 5. Stall force as a function of the frequency of the dichotomous driving and 
prediction of equation (fT0|) . The insets show, for four values of F p , n°^ ot (v) (upper 
inset) and the mean velocity (lower inset) computed (symbols) and predicted (lines) 
by equation ([Tj) . 



Assuming again a weak dependence on v and F p for n™ ot we can rewrite the previous 
equation as 

v ~ <R + <R (7) 
1 + vq/u l + v/v 

where v p LR = —F p /N is the low frequency limit of v and v p HR the high frequency one, 
v P hr = (F M n™ t;HR - F p )/N. 

Figure El shows (lower inset) our numerical results for the polymer velocity in the 
presence of a pull and it allows for a comparison against the theoretical predictions. 
As seen in the figure, the polymer mean velocity can be understood in terms of 
equation (J7J). Figure El shows the excellent agreement between this equation and 
the numerical calculations once the low and high ATP concentrations (low and high 
frequency) values are determined. The value of v p HR present in that equation can 
be evaluated experimentally by using the values of the velocity curves for high ATP 
concentrations [15J. A rough estimation can be also given by v p HR ~ (FmLm / do — F p ) /N . 

We have also studied the frequency dependence of n^ ot (^; F p ). As visible in the 
top inset of figure El this quantity depends weakly on v for all the pulling forces and 
decreases slightly with F p for all the frequencies. Then the velocity of the polymer 
shows the expected Michaelis-Menten type behaviour for all the different pull forces 
used in the calculations (F p = 0.1,0.22,0.34). We have verified numerically that we 
can consider n™ ot (v,F p ) ~ n™ ot KR (F p ) ~ n mot;HR,(0) — a ^p- I n our case a — 2-35 and 
n mot-HR.(0) — 5.59, obtained through a fit procedure of w^ ot . HR (F p ), the high rate values 

<rf F p ). 
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Stall Force. We will study now the stall force F sta u of the system, which is the value 
of the pull force F p for which the polymer velocity is zero. In formulas 

f stall — £p{V — 0) — r M — 7—. / — r M '-. : j • (o) 

1 + Vq/V 1 + Uq/U 

Here we have used the weak frequency dependence of n^ ot as suggested by the inset of 
figure El A simple approximation for n^ ot .- HR (F p ) is to assume n^ ot . HR = Lm/^o- As 
discussed above, a more realistic one is to take the linear dependence n™ ot HR (F p ) ~ 
n mot;HR(0) — aF p . Substituting the latter expression in equation (jHJ), we obtain 

„ _ p n mot;HR(°) „ f Q \ 

1 + vq/v 1 + vo/v 

By solving with respect to F p , and introducing F sta ii, we obtain 



rpHR 

1 stall 



(10) 



1 + 6i/ / ^ 

where b = 1/(1 + aFjy) and F^ h ~ 6FMn^f ot . HR (0) is the high frequency stall force. 

Figure [5] shows our numerical results for the polymer stall force problem and 
compares them to our theoretical predictions. We find an excellent agreement with 
the predictions of our equation ([ID]) , where b ~ 0.81 and Fj^ ~ 0.453. These numbers 
results from the values of a ~ 2.35 and ^mot-Hii(O) — 5-59 obtained above.. 

In this model, the stall force is not weakly dependent on the frequency u, but 
changes from to 4.5Fm in the range of rate variation of figure [5j Thus an experimental 
investigation can easily verify the stall force frequency dependence by lowering the ATP 
concentration in the surroundings of the motor. This way the working model can be 
verified and tested, and also compared with different models, [j] 

Stall force for different k. We have made a number of simulations to study the 
dependence of the stall force with the elasticity constant k. The results are plotted 
in figure [6] where the stall force is drawn as a function of the frequency v for various 
values of k. We can observe that the stall force increases with k and shows a saturation 
behaviour at high k (see the curves for k = 5 and k = 10) similar to the one of the 
velocity (inset of figure [3]). 

The MM behaviour predicted by equation ffTUl) is observed for all the curves. The 
inset of the figure shows the behaviour of the coefficients F s ^ l and b of the equation. 
As expected, the rigid chain (high k) follows an exact MM law, b = 1 in equation (ITOj) . 

Polymer packing. One of the recent most relevant activity in the field is the study 
of the translocation features in the DNA packing problem driven by molecular motors 
[T2| [14j [15]. One example is the 029, a bacteriophage virus which is able to inject his 
DNA in a bacteria in order to replicate, and then repack it in his capsid. Remarkable 
experiments [H] have measured the force of the motor as a function of the number 

| We refer in particular to the outcomes of a similar motor model assisted by both a sinusoidal time 
dependence [TB], and a pure dichotomous driving [T7]. Besides the not trivial behaviours of Fstaii as a 
function of the frequency, a very weak variation of its values is there observed. 
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Figure 6. Stall force F^ ta[l for various values of the elastic parameter k. Lines are the 
best fit according to equation (|TU1) and the inset shows the values of the fit parameters 
for the different k. 



of monomers entered in the capsid. The model here depicted is able to qualitatively 
reproduce the results there reported. 

We performed a set of calculations where the chain is pushed into a limited portion 
of space. Figure [7] shows the scheme of the motor of length Lm that pushes the polymer 
inside a region of length Leaps- The polymer starts with the last monomer on the right at 
the position x = 0. This particle (iV-th monomer in the figure) cannot pass the wall on 
the right. Once there, the polymer begins to compress and the pressure on the particles 
inside the capsid increases more and more avoiding in some cases the completion of the 
translocation process. To consider excluded volume effects and maintain the relative 
order of the monomers, a repulsive only Lennard- Jones potential has been taken into 
account in the model: Vu(r) = 4e[(^) 12 — (^) 6 ] for r < 2 1 / 6 a, and otherwise. We used 
e=l, and cr = 0.1. 

Figure [S] plots F in , the mean force acting on the last monomer entering in the 
capsid, as a function of the total number of monomers entered. The capsid width is 
Leaps/do — 8 and the chain length is iV = 36 monomers. It shows that, until the size 
of the entered polymer approximately equals the size of the capsid, the force grows up 
very slowly. Once the last monomer (JV) touches the wall of the capsid, the force inside 
grows rapidly with a saturating trend at high number of monomers. The inset c) of 
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Figure 7. Scheme of the motor pushing the polymer inside the capsid of length Lc, 




Figure 8. Force on the polymer as a function of the number of monomers entered 
inside the capsid (iVj n ). This curve resembles the experimental outcome reported in 
|14j . The length of the capsid is Leaps = 8 and the number of monomer is N = 36. 
Insets: a) trajectory of the centre of with b) a small portion showing stops; c) 

chain velocity as a function of iVj n . The activation rate of the motor for the calculations 
presented in the insets is v = 1. 



the figure plots the polymer mean velocity vs the number of monomers in the capsid. 
We find a behaviour similar to that experimentally observed in [14J. It shows that the 
model here introduced can depict qualitatively the packing features of the 029 motor 
despite of its simplicity. 

Biological values. The model we present here is a simplification of both a real 
polymeric chain and a molecular motor: in the example here used the DNA and the 
motor of the 029 bacteriophage. 

The model is in principle adaptable to any translocation process mediated by 
ATP-based motors, provided that a proper scaling of the measurable quantities can 
be done. Experimental possible value of the working time of the motor is Tm = 10ms, 
of the affinity constant hu ^ 30/iM [12], and of the saturation velocity v max ^ 103 
DNA base pairs/s = 70nm/s [15]. The saturation velocity sets the vhr value in 
our model, and the frequency at the Michaelis-Menten concentration corresponds to 
v = i/q = 1 /Tm- Finally, the stall force value found experimentally in the high frequency 
limit is Fg+QU = 57pN [13] . Using this value we can estimate Fm ~ 12pN. In the same 
way the energy consumed per cycle can be obtained for large ATP concentration as 
W = n™ ot v max F M T M . With the parameters used, W ~ 12kT room , which is less than the 
energy provided by an ATP molecule, i.e. around 20kT room . 
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5. Discussion 

In this work we have introduced a simple model for a molecular motor which consumes 
ATP and produces mechanical work pushing a polymer chain. 

Other simple models have been used to describe translocation features of DNA. 
For instance, in the works by Downton a Id polymer chain jointed together by a FENE 
potential is studied [HI [19]. The chain moves aided by a flashing extended ratchet 
potential. By contrast, our potential acts in a limited region of space which is a more 
realistic approach to translocation through a pore. 

In our work, we calculate the mean translocation time and the velocity of the 
polymer as function of the activation frequency of the motor. The latter shows a very 
good agreement with the experiments and follows a clear Michaelis-Menten dependence 
with the ATP concentration. We see that such MM laws arise in a natural way in the 
description of the system as a consequence of the kinetics of a machine which remains 
active for a certain given time. It appears that the working time average of the molecular 
motors which use ATP, is actually the origin of the Michaelis-Menten law in ATP-motor 
assisted dynamics. 

We have also studied the behaviour of the polymer in the presence of a pulling force 
and obtained analytical expressions for the stall force of the system as a function of the 
ATP concentration. Such expressions, a corrected MM equation, show an excellent 
agreement to our computed results. Finally, the force inside the capsid and the velocity 
of the chain as a function of the amount of polymer packed have been also evaluated 
finding good qualitative agreement with the experiments. 

We have studied a one dimensional model. Polymer translocation experiments are 
usually performed with the help of optical traps. That way, the polymer (DNA or 
RNA) is held almost completely stretched out. With respect to the capsid effect, the 
confinement introduced in the model describes a saturating pressure inside the capsid 
without taking into account specific geometrical details and polymer recoil effects. For 
these two reasons the one-dimensional model is a good enough first approach to study 
the process. 

Besides the fact that technological improvements allow nowadays very precise 
investigations at the nanoscale, the microscopic working details of the bacteriophage 
motors are not completely understood. The model is close to one of the two mechanisms 
recently proposed in [12], where the force acts by means of steric interactions, without 
chemical bonds with the DNA. Since DNA packing motor 029 is a well studied motor 
protein with several intriguing and unexplained features we have compared our results 
to experiments with this motor. In general we have found a good agreement. However, 
a detail model for 029 should include many other aspects and it is beyond the spirit and 
purpose of our work. Between these aspects we could mention a precise modelling of 
the inactive stage, the inclusion of more complex waiting time statistics or the addition 
of existing structural data. 

The results here presented describe general features of a motor which actuates with 
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a mean force during its cycle. In this sense the qualitative results here reported need 
some adjustments of the parameters if applied in concrete cases and can be applied to 
a wide class of ATP-based motor, independently on the inner functioning of the very 
motor. Our work shows that our model could be a good starting point to develop 
detailed descriptions of different motors. 
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